
*This is the replication file for the models in: "Singh, Shane P., and Ryan E. Carlin. 2015. "Happy Medium, Happy Citizens: Presidential Power and Democratic Regime Support." Political Research Quarterly 68(1): 3-17.

*Model 1
xtlogit dembest c.negrettoilp##c.negrettoilp lngdppercappppWB gdppercapgrowthWB durable age educ female egoretro polint proximity crimevictim , i(cntryyearnum) intmethod(ghermite)
*Model 2
xtlogit dembest c.polcon3##c.polcon3 lngdppercappppWB gdppercapgrowthWB durable  age educ female egoretro polint proximity crimevictim , i(cntryyearnum) intmethod(ghermite)
*Model 3
xtlogit dembest c.execbox2##c.execbox2 lngdppercappppWB gdppercapgrowthWB durable  age educ female egoretro polint proximity crimevictim , i(cntryyearnum)   //* did not use Ghermite method for Model 3 because Stata drops GDP Growth and fails to converge when we do.

**for Models 4-6, you must install "gllamm" ("xtologit/xtoprobit" and "meologit/meoprobit" were not yet available when we wrote the paper). // Here we use "hand-made" squared versions of the executive powers measures, as gllamm doesn't recognize the "##" syntax.

*Model 4
gllamm demsat negrettoilp negrettoilpsq lngdppercappppWB gdppercapgrowthWB durable age educ female egoretro polint proximity crimevictim , i(cntryyearnum) link(ologit) family(binomial) 
*Model 5
gllamm demsat polcon3 polcon3sq lngdppercappppWB gdppercapgrowthWB durable age educ female egoretro polint proximity crimevictim , i(cntryyearnum) link(ologit) family(binomial) 
*Model 6
gllamm demsat execbox2 execbox2sq  lngdppercappppWB gdppercapgrowthWB durable age educ female egoretro polint proximity crimevictim , i(cntryyearnum) link(ologit) family(binomial)

